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Abstract 

The one-dimensional harmonic oscillator wave functions are solutions to a Sturm- 
Liouville problem posed on the whole real line. This problem generates the Hermite 
polynomials. However, no other set of orthogonal polynomials can be obtained from 
a Sturm-Liouville problem on the whole real line. In this paper we show how to char- 
acterize an arbitrary set of polynomials orthogonal on (— oo, oo) in terms of a system 
of integro-differential equations of Hartree-Fock type. This system replaces and gen- 
eralizes the linear differential equation associated with a Sturm-Liouville problem. 
We demonstrate our results for the special case of Hahn-Meixner polynomials. 



*PACS numbers: 02.60.Nm, 05.30.Fk, 21.60.Jz 



1 Introduction 

In this paper we will examine the system of N coupled nonlinear integro-differential equa- 
tions defined on the whole real line 



w-i^' + ^w + "- w 



1 N-l ^oo 

«»W-;E / dyu*(y) 

Z ■ n J — OO 



j=o J -°° 



TT'(x) - U'(v) 

x v ' _ yyj [u i {x)u j {y)-u j {x)u i {y)} =0 (z = 0, 1, 2, . . ., N- 1). (1) 

x y 

This system is an eigenvalue problem in which there are N eigenfunctions U{(x) and N 
corresponding eigenvalues /j,\ • The eigenfunctions are constrained by the orthonormality 
condition 

/OO 
dxui(x)uj(x) = 5ij. (2) 

-OO 

The single function U(x) together with the integer N completely characterize the eigenvalue 
problem under consideration. The function U(x) is restricted by the requirements that it 
be twice continuously different iable, that it be bounded from below, and that it satisfy the 
inequality constraint 

/OO 
dxe~ u{x) <oo. (3) 

-OO 

This eigenvalue problem constitutes a system of Hartree-Fock equations that describe 
the ground state of a one-dimensional Fermi gas having iV particles. The functions Ui{x) 
are the single-particle states of the Fermi gas and the numbers fx\ are the corresponding 
chemical potentials associated with (|2]). 

As will be shown later, (|1]-[3|) arise in the study of random Hermitian N- dimensional 
matrices $ whose statistical weight is given by [I|] 

^($) = I e -^(*)rt, (4) 

Zj 

where Z is a normalization factor |J . Associated with such an ensemble of random matrices 
is a set of polynomials Vi(x) that are orthonormal on the whole real line with respect to 



the weight in (§ 0: 

/oo 
dxV i (x)V j (x)e- u(x) = S i3 . (5) 

The relation between the polynomials {Vi(x)} and the eigenf unctions {ui(x)} of (|l|) is 
simply 

Ui(x) = e-* u MPi(x). (6) 

If we substitute (|j) into ([!]) and use fl2|) we can calculate the eigenvalues {/4 }. The connec- 
tion between the eigenvalue problem fll]) and the theory of random matrices and orthogonal 
polynomials will be discussed in more detail in Sec. III. However, the principal consequence 
of this connection is that (|j) is the unique solution [u (x), u±(x), u 2 (x), . . . ,Un-i(x)] to 
([I]) and (|2|) for a given function U(x). Physically, this uniqueness is a consequence of the 
Fermi gas being in its ground state. To the best of our knowledge, the eigenvalue problem 
in (Q) is new in the general theory of orthogonal polynomials. 

In the special case U(x) = \x 2 , (0) simplifies to a set of N decoupled linear differential 
equations: 

/ ' ]i +-x* + 2^-N+)\u l {x) = V. (7) 



V dx 2 4 n 2 / 

These are the Schrodinger equations for the lowest iV eigenstates of the one-dimensional 
harmonic oscillator from which we immediately infer that /4 = |(iV — 1 — i). In the case 
of the harmonic oscillator, /4 is the depth of the ith state below the Fermi level. 

We will now show that the Hermite polynomials are the only set of orthogonal polyno- 
mials that are determined by a Sturm-Liouville problem on the whole real line. Let {Vi(x)} 
be a set of orthogonal polynomials satisfying (||). Assume further that the functions {ui(x)} 
as defined in @ are the eigenfunctions of the general Sturm-Liouville problem Q 



_d_ 
dx 



p(x)—Ui{x) 

dx 



+ [q(x) — HiW (x)]iii(x) = (— oo < x < oo), (8) 



where p(x) > and q(x) are continuous functions and {//«} are the associated eigenvalues. 
Since the domain of (|8[) is infinite, the only boundary condition imposed on the eigen- 
functions {ui(x)} is the normalization condition @. We emphasize that p(x) is strictly 
positive for all x on the whole line.Q Moreover, since the Sturm-Liouville problem in (||) 
is one dimensional, the orthonormality conditions @ imply that the spectrum {/Xj} is dis- 
crete, nondegenerate, and bounded from below. Following (0) we note further that the 
weight in (H) is fixed by W(x) = 1. Substituting the functions {ui(x)} from @ into the 
eigenvalue problem (ffl), we obtain 



p(x)V?(x) + 



w'(x) . . ,. . 

2—^-p(x)+p'(x) 



V'(x) + 



, ,w"(x) ,, s w'(x) , N 

p[x) — -— + p \x) — — + Hi — q{x) 



w(x) 



w(x) 



Vi(x) 



= 0, 

(9) 



where w(x) = e~^ u ^ x \ Equation (Q) holds for all real x and may be used to express 
p(x), q(x), w(x), and the spectrum {/ii}i>3 in terms of the first three polynomials and 
their corresponding eigenvalues as we now show. 

Let the first three orthonormal polynomials in M) be 



V (x) 



1, Vi(x)=ax + b (a^O), V 2 (x) = ex 2 + dx + e (c ^ 0), (10) 



where we have assumed the normalization 



dx w 2 (x) = 1. 



Substituting P an d V\ from flip]) into (|j) we obtain 

, . , ,w"(x) ., ,w'(x) 
q(x) = p{x) — -— + p [x)— rT - + fiQ, 



w\x) 



n w'(x) , s ,, > 
2— ^o(a;) +i/(x) 



W X 



10 (x) 
+ (jii-f*o)lx+-\=0. 



(11] 



(12) 



-'^In our case the domain of definition of the Sturm-Liouville problem is the whole real line, — oo < x < oo. 
In cases where the domain of definition is a compact subset of the real line, p(x) may vanish at its boundary 
points, which then become singular points of the Sturm-Liouville problem. In the case considered here the 
function p(x) may vanish only at x = ±oo. 



These results allow us to rewrite (If) as 



p{x) V'lix) - (jn-fio) I X + -J 7*(z) + (ja-fJto)Vi(x) = 0. 



Substituting Vzix) from (|10|) into (y|) we find that 

p{x) = 



H2~ Ho \ 2 , 
Hi~ Ho o — I x 



0*i - no) - - 0*2 - A^i) ^ 



:r 



+ 



2r: 



bd 



(A*i - A*o) (A*2 - A*o) e 

a 



(13) 



(14) 



The function p(x) is strictly positive for any real value of x. Thus, either p(x) is a 
quadratic polynomial and the requirement of positivity implies the two inequalities 

A*2 + Ho - 2 A i i < 0) 



,2 / b d 



{Hi~Ho) 7T + (H2 ~ Ho) (H2 + Ho ~ ^Hh 

\a 2c 



£ 

2c. 



< 0, 



or else p(x) is a positive constant and we have 



H2 ~ Ho 







Hi Ho 
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(Hi 


~ Ho) (A*2 - 
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. d 
-^2c 


1 
2c 


0*i- 


\ bd , 

-Ho) (H2 - 

a 


- Ho) a 



0. 



0. 



> 0. 



(15) 



(16) 



Equation ( |14D is only one of many possible expressions for p(x). Indeed, one may 
alternatively solve ([H|) for p(x) in terms of any one of the polynomials {"Pj}j>3 as 

(x + J) (/ii - /i ) 7^(x) - {Hi ~ Ho) Vi{x) 



p(x) 



V'({x) 



(17) 



which owing to ([14]) , must be either a quadratic polynomial in x or a positive constant 
for all i. This is precisely the condition that determines the whole spectrum in terms of 



fjLo, fjLi, and fj,2 • To see this, we take the limit of (|TT| ) as |x| — >■ oo and using (|T3]) we find 
that 



/'; 



M2 -yWo 



(/*i - ^o) 



* 2 + 



2 0*i - Mo) - 



A*2-A*o 



Z + /X0- 



is) 



The condition that the spectrum in ([18]) increases monotonically implies the inequality 



/i 2 + no- 2/ii > 0. 



(19) 



If p{x) is a quadratic polynomial, then the inequality /x 2 + Mo — 2mi > must hold, but 
this contradicts the first equation in flT5|). It also contradicts the second inequality in flTSP 
because if both these relations are valid, they imply the inequality 

2 



2c, 



<0. 



(20) 



However, from ([K]) this is precisely the condition that Vi{%) be nonvanishing on the real 
line, which implies that its sign never alternates along the real line. This conclusion con- 
tradicts the orthogonality of Vi and Po-Q Therefore, (|T5|) can never hold, and (0), which 
implies that p(x) is a positive constant p, is the only possibility. Finally, using (|TT|-|r2|) we 



find that 



and 



w 2 {x) = e 



-U(x) 



2np 



1/2 



Ml — /Xq 



e 2 p 



/ x (wi-Wo) 2 / . b\ 3 1 



(21) 



(22) 



for which (H) becomes the Schrodinger equation for an harmonic oscillator whose equilib- 
rium point is at xq = —b/a, and the polynomials {Vi(x)} are the corresponding (shifted) 

2 There is another way to see that (|20| ) is false. We note that the coefficient e in V2{x) in ( jlOl ) is actually 
not arbitrary; it may be determined in terms of the coefficients a, b, c, and d as follows. We consider the 
following three orthonormality relations: (1) f dxw 2 (x)T'o(x)'Pi(x) = 0; (2) J dxw 2 (x)'PQ(x)T > 2{x) — 0; (3) 
J dx w 2 (x)Vi (x) — 1 . These equations may be regarded as three linear simultaneous equations for three 
unknowns, the first and second moments of w 2 (x) [the zeroth moment is given in (|llj)] and the coefficient 
e. The solution for e is e = (abd — c — b 2 c)/a 2 . When this value is substituted into the inequality in (|20|), 
the inequality takes the form [(ad — 26c) 2 + 4c 2 ]/(4a 2 c 2 ) < 0, which is manifestly impossible. 



Hermite polynomials and the eigenvalues are //, = (/ii — /j,o)i + n$. This concludes the 
demonstration that the Hermite polynomials are the only set of orthogonal polynomials 
that are determined by a Sturm-Liouville problem on the whole real line. 

This paper is organized as follows. In Sec. II we present a proof of the integro-differential 
eigenvalue equation (fj) based on a variational calculation. In Sec. Ill we give a general 
formula for the eigenvalues associated with an arbitrary system of polynomials having 
an even weight function. This formula for the eigenvalues is useful because it gives the 
spectrum for a one-dimensional system of interacting fermions. Finally, in Sec. IV we 
illustrate equations ([[]) and (Q) for the special case of Hahn-Meixner polynomials. We 
calculate the first few eigenvalues and find some of their general features. 



2 Orthonormal Polynomials on the Real Line 

In this section we use a variational principle to prove that polynomials orthonormal relative 
to an arbitrary weight w 2 (x) = e~ u ^ x ' on the whole real line must obey the system of 
Hartree-Fock equations in ([!]) J2]. We show that as a consequence the functions in (|(|) 
are the unique solution of (H). [Recall from Sec. I that the system of equations in (|l|) 
reduces to the trivial set of uncoupled Schrodinger equations (|7|) for the eigenstates of an 
harmonic oscillator (i.e., the case of Hermite polynomials) when the weight is Gaussian 
w 2 (x) = e~ x Z 2 .] 

Consider a quantum mechanical system whose degrees of freedom are elements of an 
N x N Hermitian matrix $. The Hamiltonian for this system is the positive semi-definite 
operator given by 



H=ilV 



4 + ^<*>)GI + ^w 



(23) 



We will require that U(Q) be a matrix potential function satisfying the restriction in (|3|). 
This clearly implies that the potential in (p3|) is bounded from below and grows to plus 
infinity as the matrix eigenvalues become infinite. Hence, the Schrodinger operator 7i in 
(|23D has a well-defined spectrum. This Hamiltonian is symmetric under the adjoint U(N) 
transformation $ — > f2$f2t Therefore, Ti possesses a unique (U(N) singlet) normalizable 
ground-state vector ^o^)- This ground state has energy zero and is given byQ 

*(>(*) = -^exp[-iTr[/($)] 



Z = f d N2 $ex-p[-Tr [/($)]. (24) 



As is well known, the Laplacian over Hermitian matrices acquires the form ||, [|, 18 



Q2 1 N d 2 

Tr— — ■ = — y^ t-jA(i) + [U(iV) angular momentum terms], (25) 

o9 2 A{x) ~[ oxj 



3 It is easy to verify that this state has zero energy by substituting (M) into (p3). Thus, since TC is 
non-negative, ^o hi (p4T) is the ground state because the latter is unique. 



where X{ are the matrix eigenvalues and A(xi) is the Vandermonde determinant. 

For a generic potential £/"($), the Hamiltonian in ( |23| ) contains long-range two-body 
interaction terms between the Fermions^ 

1. 



T^int = — jTr 



<l) c, ' ( * ) 



1 v gfe) - U'ixj) 

A ^ r- — t ■ 

^ i,j a* X J 



Thus, a one-dimensional collection of eigenvalues may be considered as an interacting Fermi 
gas. 

The eigenstates of 7i have definite U(N) quantum numbers and thus fall into definite 
U(N) representations. Because of this symmetry it is useful to introduce matrix polar 
coordinates $ = QXQ^ , where fl G U(N) and X is a real diagonal matrix whose entries, 
Xi, X2 ■ ■ ■ ,xn, are the eigenvalues of $. The transformation from Cartesian matrix coor- 
dinates <frjj to polar matrix coordinates fl^, and Xij = Xi5ij is associated with a Jacobian, 
which is proportional to A 2 (xi). Under this transformation 7i and \l/o transform as 

H -> AHA' 1 , 
^o -> A^o- (27) 

This implies that the transformed self-adjoint Hamiltonian is symmetric and that the trans- 
formed ground state is totally anti-symmetric under the interchange of any two of the 
eigenvalues Xj. These symmetry properties remain valid for any U(N) singlet eigenstate of 



7i. Therefore, the singlet sector of ( p3]) is equivalent to a one-dimensional Fermi gas with a 
fixed number N of particles [IS | , where the eigenvalues xi are considered as the coordinates 
of the N Fermi particles. 

Let us consider the class of singlet states of the form 

V v {$) = exp[-^FrV{$)], (28) 



4 Note that if U is a polynomial of degree less than or equal to three, there are no two-body interactions 
in @. 

9 



where tyy is assumed to be normalizable. In matrix polar coordinates this becomes 



*y(x<) = Afdetij (x) x ) 



exp 



1 N 

fc=i 



(l<i,J<iV), (29) 



where we have used the identity A(xj) = detjj (x* -1 ) and J\f is a normalization constant. 
It now becomes clear how this discussion of a Fermi gas relates to the theory of poly- 
nomials H]. First, note that 

det.j (x)- 1 ) = det^ [Pi-ifo)] (1 < i,i < iV) (30) 

for any set of monic polynomials Pi(x) = x l + (lower powers of x). Thus, the wave function 
(J29"D is a Slater determinant [p^j : 

^y(^) = -^=detij [vi_i(xj-)] , (31) 

where Vi(x) = v i P i (x)e~ v( - x ^ 2 are interpreted as normalized single particle states in the N 
fermion wave function fl3T|), which will be used as a trial wave function in a variational 
(Hartree-Fock) calculation of the ground state of the Fermi gas and M is a normalization 
constant. 

Therefore, these functions must satisfy the orthonormality conditional 

dxe- v{x) v i (x)v j (x)=5 ij} (32) 

which implies that the polynomials {Pj(x)} in (^) are the set of polynomials orthogonal 
on the real line with respect to the weight e~ v( - x \ In particular, for V(x) = U(x), we have 
^v( x ) — ^u( x ) — ^o( x ) an d {vi(x)} are identical to the orthonormal functions {ui(x)} 
defined in (|6|). 

Consider the expectation value of the Hamiltonian Ti in fl23|) in the normalizable state 
of©: 

1 r ]N 2 ^-hTTVnj^-hTrV 



E (V) = ■£- / d" Qe-^-He-? 1 ™, (33) 

Zy J 

10 



where 



Clearly, 



Z v = / d N2 $e- TrV = {V V \V V ) < oo. 



S {V) > 



(34) 



(35) 



for any normalizable \l/y. 
Since we have 



£o(U) = 0, Z v < oo, 



(36) 



we expect ^u = ^o to be an absolute quadratic minimum for £q(V) in the space of functions 
of the form given by (p8|) and Q34D . Indeed, for V infmitesimally different from [/, 



V = U + 8U, 



(37) 



we find that 



So(U + 5U) = {V v \-Tr{6U'Y\® v )/Z v + 0{5U) 6 > 0. 



(3* 



Expressing ( |38|) in matrix polar coordinates and using the Slater determinant representation 
( |3~1~1) for \&y (and for \&i/), we find, after some trivial manipulations, that 



iV-l „ 

£ (U + SU) = Y, dxv* 

i=0 



X 



w-i^» + ^w + "- w 



^i(x) 



-7 E iJ / dxdyv*(x)v*(y) — — —[v^v^-v^v^y)} - ^ ^f ] 

4 J x v «=o 



1 N ~ 1 r 

8 i=0 J 



(39) 



Here, the {/4 jilo 1 are a se ^ of iV Lagrange multipliers (also known as chemical potentials) 

that will enforce the unit normalization condition on the {vi} in the variational calculation.^ 

5 It is unnecessary to include in (p9) Lagrange multipliers to enforce all orthonormality conditions (|32| ) 
because the {v{\ which result from the variational calculation turn out a posteriori to be orthonormal p9[ . 
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Taking the variational derivative of fl39| ) with respect to v*(x), at V = U, we find that 
SSo(U + SV)/5v*(x)\ su=0 = {[-\dl - -U" + V + f4*]v t (x) 



--EJ / dyv*(y) [vi(x)vj(y) - v j (x)v i (y)]}\ Vi=Ui 

2 J J J x — y 

= 4E J dy\uM\ 2 SU\y)y^\su=o = 0. (40) 

This expression must vanish because 



r _, . 5£/'(z) 
ovi{y) 

must be a regular kernel in function space for V infinitesimally close to U because V{x) 

defines the Vi(x) uniquely, and vice- versa [V(x) = —2 log (^p )]■ At the point V(x) = U(x) 

in function space the set of functions {vi (x)} coincides with the set of functions {ui (x)} 

defined in (H), where the {Vi(x)} are the set of polynomials orthonormal with respect 

to the weight e~ u ^ on the real line. Moreover, uniqueness of the ground state of the 

quantum-mechanical system defined by ( p3| ) implies that (||) is the unique solution to (|1|). 

This completes our proof of the assertion made in Eqs. (|1|) and (||). 

We conclude this section by proving that the first v orthonormal polynomials, Vq, Vi, 

Vi-, ■ ■ -i Vvi associated with the weight w 2 (x) = e~ u ^ x ' satisfy the useful sum-rule 

u ~ 1 ™ '- U'(x)-U'(y), 



r°° r°° I J It) — I J I11) 
£ / dx w\x) / dy w 2 (y)^ ^Pn(x)P m (y) [V n {x)V m {y) - V m {x)V n {y)) 

. r, J —CO J — CO X y 

"- 1 roo 

2J2 dxw 2 (x)[P' n (x)] 2 . (42) 



n=0 " ~°" 

To prove (^) we let $ be a v x v Hermitian matrix and define 

/^ 2 $e- TrC/ (*)Tr[(^) U'($) 



S„ 



Integrating by parts over $ we obtain straight forwardly 



(43) 



Jcf 2 $e- T ^>Tr[[/' 2 ($)] 
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In terms of the Slater determinant ( |3T| ) (with V = U) the latter equation becomes |T9[ 

v- 1 /-oo 

S U = J2 dx w 2 (x)V 2 n (x)U' 2 (x). (45) 



n=0" 



On the other hand, using the identity 



Tr 



m) u ' w 



" U'(x m ) - U'(x. 



E 



m,n=l 



Xm X n 



-=L (/ ( x ™) + E — i — i ( 46 ) 



m=l 



m^n=l 



*^m *^n 



and the same Slater determinant as above, we obtain [19| an alternative expression for S v 



as 



s v 



v-\ 



oo 



, ■- r°° U'(r) — U'(v) 
£ / dxw\x) dyw 2 (y)^ m V n {x)P m (y) \V n {x)V m {y) - V m {x)V n {y)\ 

m ,n=o J -°° ■ y -° o x ~y 



oo 



/oo 
cW(x)7^(x)£/"(x). 



(47) 

Equating © and (0), using the identity e~ u ^[U' 2 {x) - U"(x)\ = £ze~ U{x) and the 
orthonormality condition fl5|), we obtain (f|I|) as required. 
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3 Eigenvalues of the Integro-Differential Equation 

In the previous section we have proved that the set of functions in (|6|) is the unique solution 
of the system of integro-differential equations in (p. The eigenvalues {/4 } of ([[]) played 
absolutely no role in that proof because we already knew the exact form fl24]) of the Fermi 
gas ground state and used it to show that ([!]) and (H) hold. However, these eigenvalues are 
an indispensable ingredient of the system (P because they are Lagrange multipliers that 
enforce the orthonormality conditions (|32|). They therefore encode important information 
about the orthonormal polynomials (|6]) in much the same way such information is encoded 
by the set of moments of the weight e~ u ^ or by the set of coefficients of the generic three- 
term recursion relations among the polynomials Vi [^, |21], [1], [3]]|] Since the polynomials 
Vi are known once the weight e~ u ^ is given, we can use (Q) to determine the eigenvalues 

Multiplying (|IJ) by u*(x) from (|B|) and integrating over x we find, using (|2|), that the 
eigenvalues are given by 

1 /■« , „ , I"" 1 



rfV = -- J^ dx w\x) {V' n f + - Y,K nm (0 < n < A - 1). (48) 

Here, as before, w 2 (x) = e~ u ^ x > and 



m=0 



K = K 



oo 



.2,^ / J ....^..^'{x)-W{y) 



dx w (x) / dy w (y)- 



oo J —oo x — y 

xV n (x) V m (y) [V n (x) V m (y) - V n (y) P m (x)] 



oo roo 



oo J — oo 



dx / dy)C nm (x,y), (49) 



where we define 

y- i \ 2 f \ 2/x/o a ,'P n (x)V m (y)[V n (x)V m (y)-V n (y)V m (x)} 

>Cnm{x, y)=w [x)w {y) (d x - d y ) . (50) 

x-y 



6 Note, however, that the latter two sets of coefficients involve all polynomials, while the eigenvalues 
{n\ } do depend on the number N of polynomials used. 
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It is obvious from (p3) that the symmetric matrix K nm has vanishing diagonal elements 

K nn = 0. (51) 

The sum-rule (^) we have found in the previous section now becomes 

N-l N-l 

E K nm = 2 E / dx w 2 (x) (V' n f . (52) 

n, m=0 n=0 °° 

Using (^)-(0) and the symmetry of the matrix K we obtain after some algebra yet another 

sum rule 

JV-l N-i roo N-l N-l i N-l 

E^-=E/ dxw\x){V' k f- E £*«+o£*u (n<N-2). (53) 

m=0 fc=n °° fe=n+l Z=0 fc,i=n 

From (|E) and (|53| ) we infer the following recursion relation satisfied by the eigenvalues 

JV-l -i JV-l 

A*?° = - E ^+ E^« (n<JV-2). 

k=n+l k,l=n 

The symmetry of the matrix K allows us to restrict the double sum on the right side of 
the last equation to the upper triangular part of K; the recursion relation then becomes 

JV-l 1 N-lN-k-l 

A#° = - E ^ N) + - 2 T, E K k , k+P (n<N-2). (54) 

k=n+l k=n p=l 

This recursion relation may be used to calculate /i^ for descending values of n holding N 
fixed. Similarly, ([P|) implies that along the N index direction there is a recursion relation 
of the form 

/^ +1) = A*? + \Ku,n (0<n<N-l). (55) 

The equations ( E5| ) and (E||) are not independent and thus they cannot be combined to give a 
single self-contained recursion relation in terms of the eigenvalues jjl)^' alone. Nevertheless, 
we can calculate p^> iteratively using the above formulas. 

The initial condition for both ([54]) and (|55|), namely the value of /%-ii m ust be calcu- 
lated separately. We show below that 

/$2i = (56) 
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for any even weight w 2 (x) = e u ( x >. 

To prove this identity we note from (pE8| ) that 

JV-l 



1 r°° 2 1 

!i = "2 j_ dx w "( x ) ( V 'n-i) + 2 E ^-i,m- ( 57 ) 



(AT) 

m=0 



Rewriting the sum in the last equation as 

N-l I N-l A7-2 \ N-l 

53 KN-l,m = I 5Z _ X! ) -^™ m ~~ X! Kmjf-l, (58) 

m=0 \n,m=0 n,m=0/ m=0 

and using the sum-rule (|5^) and the symmetry of the matrix K, we find that 

N-l roo 2 N-l 

53 K N -i, m = 2 dx w 2 {x) (V'n.i) - £ KN-l,m- (59) 

m=0 °° m=0 

Therefore, 

N— 1 <.oo . . 2 

2 #AT-l,m = / dx W 2 (X) [V'nA (60) 

m=0 • / -°° 

which upon substitution into (|57D, yields the result fl56|). 

We now give explicit general formulas for the first few eigenvalues fij^' . Our procedure 
is to use (f49|) and (|50|) to obtain a sequence of formulas for K„ m and thus for the eigenvalues. 
These formulas are obtained from the polynomial recursion relation |], [3], |2(| pi 



^n(^) = yRn+l V n +l(x) + yi2„ TV-lfa), (61) 

where the recursion coefficients R n are uniquely defined by the weight function w 2 (x). Note 
that this recursion relation generates normalized polynomials V n (x)'- 

V (x) = 1, 

X 



VRi 
x 2 -R x 



Vi{x) 

V 2 (x) 

V 3 (x) = rrr~rr^ 



V-R1-R2 

x 3 - (gi + R 2 )x 
V-R1-R2-R3 
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„, . x 4 -(R 1 + R 2 + R 3 )x 2 + R l R 3 

Pi[x) = = , 

V-R1-R2-R3-R4 



PM 



x 5 - (R 1 + R2 + R3 + R^x 3 + (R1R3 + RiR± + R 2 R*)x 



V 7 R1R2R3R4R', 



o 



V 6 {x) = j==^== =[x 6 - (R, + R 2 + R 3 + R 4 + R 5 )x* 

V -T11-K2-K3-K4-K5-K6 

+(i?ii?3 + -R1-R4 + -R1-R5 + R2R4 + R2R5 + RzR§) x ~~ -R1-R3-R5]; 

P 7 (x) = * =[x 7 - (ifr + i? 2 + i? 3 + i? 4 + ^5 + fle)s 5 + (i?iig 3 

V til til ti?, -^4 -K5 ti& 

+-R1-R4 + -R1-R5 + R\R% + R2R4 + R2R5 + R2R6 + R3R5 + R3R6 
-{-R^Rftjx — (R1R3R5 + RiR^Rft + R\R/±Rq + i^-R^,-^)^.!- (62) 

Inserting the above recursion relation into (0) repeatedly allows us to cancel the quan- 
tity x — ?/ in the denominator. 7 The ortho normality relation 

dx w 2 (x)V n {x)V m {x) = 5 nm (63) 

then simplifies the resulting expression dramatically. To present our final results we define 
a set of polynomials Q^'(x) that are conjugate to the polynomials V n (x): 

Qi N \x) = 1, 

JN), X 



Qr>(x) 



VR 



N-l 



^{N), s ^ — Rn-1 

Q 2 >(x) - 



>W 



QT'(x) 



\/Rn-iRn-2 

x 3 - (Rn-i + Rn-2)x 

^/Rn-iRn-2Rn-3 
X — (Rn-1 + Rn-2 + Rn-3)% + Rn-iRt 



(AT) , . X — {ItN-1 T- -rt,AT-2 "I" -TtAr.aJX i" JXjV-l-n-AT-3 

^ 4 ^ ~~ / p p p p = ' 

V tiN-ltiN-2tiN-3tiN-A 

Qf\x) = j= = J = = [x 5 - (iZjy-i + ifrv- 2 + i?iv-3 + R N -i)x 3 

V tiN-l-tiN-2-tiN-3tiN-4 : tiN-5 

+ (Rn-iRn-3 + RN-lRN-i + Rn-2Rn-i)x], 



7 For I to — n| = 1 in ((t9[) . (|50|) this iterative process yields the Darboux-Christoffel formula. When 



|to — n| > 1 one obtains generalizations thereof 20, El 
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Qe (x) = , zj ^=[ x - (Rn-i + Rn-2 + Rn-3 + Rn-4 

V KN-iUN-2-tiN-3-tiN-4-tiN-5-tiN-6 

+Rn-5)x + (Rn-iRn-3 + Rn-iRn-4 + Rn-iRn~5 + Rn~2Rn~4 

+Rn-2Rn-5 + Rn-3Rn-5)x — Rn-iRn-3Rn-5\, 



1 

C ^fo' ~ (-^iV-1 + ^iV-2 + #iV-:i " V- : 

V ttN-l-KN-2-tiN-3-tiN-4-tiN-5-tlN-6 



Q?\x) = /n = „ \ „ ==[x 7 -(i? jV _ 1 + R N - 2 + Rn-3 + R 



+Rn~5 + Rn-6)X + (Rn-iRn-3 + Rn-iRn-4 + Rn~iRn~5 
+Rn-iRn-6 + Rn-2Rn-4 + Rn-2Rn-5 + -RjV-2-RiV-6 + RnsRn-5 
+RN-3RN-6 + Rn-aRn-g)^ — (Rn-iRn-3Rn-5 
+Rn-iRn-3Rn-6 + Rn-iRn-4:Rn-6 + -^-2-^-4-^-6)^] 5 (64) 

where we have replaced i? n in (B2I) by Rn-u- For each value of iV this set of polynomials is 
orthonormal as can be seen from the fact that they satisfy a three-term recursion relation 
whose recursion coefficients are Rn-u- 
We then obtain the formula 

K n ^ k = -== / dx w 2 {x)V' n {x)QT\{x)V n - k {x). (65) 



Once we have determined the matrix i^ nm we can then use (p5|) to construct a formula 
for the eigenvalues /x^: 

i N-n-l 

/*?° = 2 E *«^-i (n<iV-l). (66) 

We have been able to find closed-form expressions for some of the K nm : 

K n , n -i = — (n > 1), 

tin 

r\ n—1 

K n , n -2 = -jt-b — E^ (n>2), 

I n-2 

^n,„- 3 = tt^ 5 — 53[i^_ 1 + 2i^ + 4i2 fc i2 fc+ i-4fl n _ii2 fc ] (n>3), 

tin-tin-l-tin-2 f.=l 



n n— 3 

Kn^n-i = 2^ [{Rn-1 + Rn-2) Rk ~ 2(-R„-l + R n ~2){Rk + 2-Rfc-Rfc+l) 

tin,tin-\tin-2tin-?, k= \ 

+3R k R k+1 R k+2 + Rl + 3R 2 k R k+1 + 3R k R 2 k+1 ] (n > 4), (67) 

and so on. These formulas become increasingly complicated as one moves further away 
from the diagonal of the matrix K. 

From (|66|) and (|6^) we obtain the following formulas for the first few eigenvalues: 

MJ2i = o (jv>i), 

*#?. = ^f 1 (N>2), 

l& = TTJT^-R 12(Ri + --- + Rn- 2 ) + (N-2)R n „ 1 ] (iV>3), (68) 

ZHn-I-tLn-2 

and so on. Again, these formulas become more complicated as the difference between the 
two indices of \x grows. 
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4 Eigenvalues of the Hahn-Meixner Polynomials 

In the previous section our discussion was completely general; it applies to the case of 
an even weight function w 2 (x) on the whole real line. To illustrate the results of the 
previous section we now focus on a special class of polynomials known as the Hahn-Meixner 
polynomials. In the context of the discussion above, these polynomials are single-particle 
states from which one constructs the ground state of a one- dimensional Fermi gas. 

The Hahn-Meixner polynomials || |7|, |8], |9], |10, [Ll] are a four-parameter family of poly- 



nomials V n (x) orthogonal on the whole real line (—00 < x < 00). These polynomials are 
expressible in terms of a generalized hypergeometric function: 

V n {x) = N n i n 3-F 2 (— n, n + a + b + c + d — 1, a — ix; a + c, a + d; 1) , (69) 

where a, b, c, and d are parameters that may take on any complex values and the normal- 
ization N n is given by 

[(2n + a + b + c + d- l)T(n + a + b + c + d- l)T(n + a + c)T(n + a + d)} 1 / 2 

(70) 



n " T{a + c)T{a + d)n\[T{n + b + c)T{n + b + d)} 1 / 2 



These polynomials are orthonormal on the real line 

/oo 
dxw 2 (x)V m (x)V n (x) =5 mn , (71) 

-00 

where the weight w 2 (x) is given by 

w 2 (x) = — r(o + ix)Y{b + ix)T(c — ix)T(d — ix). (72) 

Hahn-Meixner polynomials appear often in mathematical physics; for example, they 
play a major role in the representation theory of the Lorentz and rotation groups |lT], |12J| . 
Recently, it has been shown that a special case of the Hahn-Meixner polynomials plays a 



crucial role in discrete-time quantum mechanics |13| , the solution of operator differential 
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equations pi, the operator ordering problem in quantum mechanics ||13||, and in the theory 



of Weyl-ordered (symmetric) operators of the Heisenberg algebra [I0|. Here we restrict 



our attention to the special case of Hahn-Meixner polynomials that appear in the latter 
applications. These particular polynomials are characterized by the parameter values 

1 3 

a = c = - and b = d = - (73) 

4 4 

in ©. 

This special case of Hahn-Meixner polynomials may be obtained in a natural way from 
the Heisenberg algebra as follows. We introduce a complete operator basis having a simple 
commutator algebra. This Hermitian basis T m>n is the sum over all possible orderings of m 
operators p and n operators q. For example, 



Ti,i = pq + qp 
To,i = q, 



T2,i = p 2 q + pqp + qp 2 , 



T 2 ,2 = (p 2 q 2 + q 2 p 2 +pq 2 p + qp 2 q + qpqp + pqpq). (74) 

These symmetrized operator products T mj „ can be rewritten as Weyl-ordered sums in the 



operators p and q with the terms in the sums weighted by binomial coefficients (TR [L6 



T mn = [m+ r ±{ n \q k p m q^ k \ (75) 

2 n m\n\ ^ Q \kJ K ' 

These basis element operators T m;T1 have particularly simple commutation properties: 

[q,T m ,n] = i(m + n)T m _i jn , 
[p,T m ,n\ = -i(m + n)T m)n _!, 

m + n + 1 

{P,T m ,n}+ = m + n+ 1 ^m+l,n- (76) 
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Thus, commuting with q and p has the effect of a lowering operator and anticommuting 
with q and p has the effect of a raising operator in the appropriate index. From the four 
basic algebraic relations (J76T) one can deduce a general formula for the commutation and 
anticommutation relations between T m>n and r r)S [[L7[ . 

The connection between the totally symmetrized operators T m>n and the special case of 
the Hahn-Meixner polynomials in ( |69"D and (|73|) is 

Tn ' n = (2n-l)!!' Pn(Tl ' l) ' 

Tn ' n+k = (n + k)\2 n +^ qk,Vn ^ Tl,1 ^ + ' ^ 

We list below the properties of the special Hahn-Meixner polynomials: 
(i) The V n (x) satisfy the recurrence relation 

nV n (x) = xPn-x{x) -{n- l)V n - 2 (x). (78) 

By comparing this special case with the general recursion relation in fl6"T|) we identify 

Rn = n 2 . (79) 

(ii) The V n (x) are orthonormal on the interval (—00,00) with respect to the weight 
function 

e -U(x) = 2/ ) = 1 (80) 

v ; 2cosh(vrx/2) v ; 

The moments of this weight function are Euler numbers: 

/oo 
dxw 2 {x)x 2n = \E 2n \, (81) 

-00 

where E = 1, E 2 = -1, E 4 = 5, E 6 = -61, E 8 = 1385, E w = -50521, and so on. 

(iii) The V n (x) have a simple generating function 

^ exp(xarctant) 

h n{X) ~ VTT¥ ■ 
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(iv) The first few V n (x) are 



V (x) 
Vi(x) 

V 2 (x) 

V 3 (x) 
V 4 (x) 
V 5 (x) 
Vb(x) 
V 7 (x) 
V 8 (x) 

Using these polynomials in 
listed in Table 1. 



x. 



2 [X 



- (x 3 - 2x 



_ ( x 4 _ 5x 2 + 3 



6 
1 

24 
1 

5! 

1 / « 35 



— [ x 5 - 10a; 3 + -i , . 



23 
Y 



(> 



x 4 + 49a; 2 , , . 



45 

6! ^v.._ a .. +49ar .__ 

— fa; 7 - 28a; 5 + 154a; 3 - 132a; 
7! V 

42a; 6 + 399a; 4 - 818x 2 + 



— x 



315 



we obtain the matrix K mv , whose first few entries are 





m=0 


m=l 


m=2 


m=3 


m=4 


m=5 


m=6 


m=7 


n = 





1 


l 

2 


l 

2 


13 
36 


7 
20 


17 
60 


1727 
6300 


n = 1 


i 





1 
2 


5 

18 


11 
36 


41 
180 


209 
900 


1201 
6300 


n = 2 


i 

2 


1 
2 





1 
3 


7 
36 


67 
300 


17 
100 


1121 
6300 


n = 3 


1 
2 


5 

18 


1 

3 





1 
1 


3 

20 


53 
300 


6017 
44100 


n = 4 


13 
36 


11 
36 


7 
36 


1 

■1 





1 
5 


11 
90 


461 
3150 


n = 5 


7 
20 


41 
180 


67 
300 


3 

20 


1 

5 





1 

6 


13 
126 


n = 6 


17 
60 


209 
900 


17 
100 


53 
300 


11 
90 


1 
6 





1 

7 


n = 7 


1727 
6300 


1201 
6300 


1121 
6300 


6017 
44100 


461 
3150 


13 

126 


1 

7 






Table 1 The first few entries of the symmetric matrix K mn associated with the Hahn- 
Meixner polynomials. The entries in this matrix were computed from (|65|). 

Applying ( |66] ) to the entries in Table 1 we can compute the eigenvalues ^ associated 
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with the Hahn-Meixner polynomials. The first few eigenvalues are listed in Table 2. 



?? 



8 9 10 



iV= 1 
N = 2 
N = 3 
N = A 
N = 5 
iV = 6 
N = 7 
N = 8 
N = 9 
N = 10 
JV = 11 



1 

2 





















3 

1 


1 

1 



















1 


7 
18 


1 

6 

















85 

72 


13 
21 


19 
72 


1 

8 















61 
15 


59 
90 


169 
450 


1 

5 


1 
10 













539 
360 


463 
600 


829 
1800 


173 
600 


29 
180 


1 
12 











286 

175 


2731 
3150 


577 
1050 


3931 
11025 


41 

175 


17 
126 


1 
11 









14713 
8400 


8083 
8400 


7331 
11760 


75703 
176400 


51403 
176400 


199 
1008 


13 

112 


1 
16 







3917 
2100 


46057 
44100 


489 
700 


21607 
44100 


20009 
56700 


145 

588 


43 
252 


11 
108 


1 
18 





346753 
176400 


22067 
19600 


134737 
176400 


877307 
1587600 


214189 
529200 


9257 
317520 


22633 
105840 


541 
3600 


49 
540 


1 
20 



Table 2 Table of the eigenvalues //^ associated with the Hahn-Meixner polynomials. 
These eigenvalues are obtained from the entries in Table 1 using (|66|). 

For the case of the Hahn-Meixner polynomials the general formulas in fl5S|) become 



..(N+l) 
H'N-l 

(N+l) 



..(N+l) 



JN+l) 
PjV-4 



U (N+1) 



JN+X) 

H-N-G 

..(N+l) 
PiV-7 



1 

2N' 

5iV- 1 

6(iV-l)iV' 

19iV 2 - 28A^ + 3 
15(JV-2)(JV-1)JV' 
171iV 3 - 633A^ 2 + 551A^ - 45 
105(N - 3)(N - 2)(N - ljiV' 
1289iV 4 - 8980A^ 3 + 19351iV 2 - 13052iV + 840 
630(A^ - A)(N - 3)(N - 2)(N - 1)N ' 
16757iV 5 - 187721iV 4 + 738847iV 3 - 1201015iV 2 + 689652A^ - 37800 
6930(iV - 5)(N - 4)(iV - 3)(iV - 2)(iV - 1)N ' 

[2(63705A^ 6 - 1048800A^ 5 + 6531490A^ 4 - 19071030A^ 3 + 25902989iV 2 

-13249854iV + 623700)] 
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/[45045(iV -6)(N- 5)(N - 4)(N - 3)(N - 2)(N - 1)N], 
H^+V = [2(72199iV 7 - 1639133iV 6 + 14772758iV 5 - 67196100iV 4 + 161413785N 3 ] 
-192719779iV 2 + 90335010iV - 3783780)] 
/[45045(iV - 7)(N - Q)(N - 5)(7V - 4)(N - 3)(N - 2)(N - l)N], (82) 

and so on. 

We have discovered that these formulas simplify considerably if they are rewritten in 
the form of partial fractions: 



2 .i!!^r 1 i) 

(JV+l) 



N' 
1 



2-3!!/^_V = t? + 



2-5!!/4 W+1) 



N N-V 

3 12 23 



3 

(JV+l) 



N N-l N-T 

15 44 107 176 



2 ^" 4 " N + N-l + N-2 + N-3' 

n «.. W+D 105 276 693 1104 1689 

2 • 9!!/47_ + 5 iJ = — + — - + — - + — o + 



(N+l) 



2.nw_V = ^r + 



AT 7V-1 N-2 N-3 iV-4' 
945 2340 4773 9360 13329 19524 



.(JV+l) 



N N-l N-2 N-3 N-4 N - 5' 

10395 24780 46767 94512 138027 185772 264207 



2-13W_V = -^ + 77-^ + 77-^ + 77-^ + ^7-^ + 



2-15!!/,r 8 1) 



AT JV-1 N-2 N-3 N-4 ' N-5 N - 6 ' 
135135 313740 566235 978480 1694655 2264940 



N N-l N-2 N-3 N-4 N-5 

2944395 4098240 

+ N-Q + N-7 ' 

(83) 

Moreover, we have found that many of the coefficients can be expressed using very simple 
formulas. For example, the leftmost coefficient is a double factorial: 1 = (—1)!!, 1 = 1!!, 
3 = 3!!, 15 = 5!!, 105 = 7!!, 945 = 9!!, 10395 = 11!!, and 135135 = 13!!. Also, the rightmost 
coefficient is the sum of inverse odd integers: 
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» - »(i + M 

176 = 7!! (l + H + 7 
1689 = 9ll (l + 5 + 5 + ? + 5 

19524 = 11!! T- + - + - + - + - + — 

VI 3 5 7 9 11 

„/l 1 1 1 1 1 1 

26247 = 13 -H 1 1 1 1 1 

U35791113 

4098240 = 15!! (- + - + - + - + - + — + — + — 
U35791113 15 



Nevertheless, we have not yet discovered a general formula for the coefficients. 



J4) 
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